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Abstract 

We study the statistics of level spacing of geometric resonances in the dis- 
ordered binary networks. For a definite concentration p within the interval 
[0.2,0.7], numerical calculations indicate that the unfolded level spacing dis- 
tribution P{t) and level number variance S^(L) have the general features. It 
is also shown that the short-range fluctuation P{t) and long-range spectral 
correlation (L) lie between the profiles of the Poisson ensemble and Gaus- 
sion orthogonal ensemble (GOE). At the percolation threshold pc, crossover 
behavior of functions P{t) and is obtained, giving the finite size scaling 

of mean level spacing 6 and mean level number n, which obey the scaling 
laws, 6 = 1.032L-1-952 ^^^j ^ ^ O.gilL^-^^o. 

PACS number (s): 77.84.Lf, 42.65.-k, 02.50.-r 
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I. INTRODUCTION 



Random matrix theory (RMT) originated from dealing with the energy levels of complex 
many-body quantum systems has become an independent new statistics.^ Different from the 
standard statistics in physics, RMT focuses its attention on the general properties of a 
number of stochastic ensembles with common symmetry. The correlation and fluctuation 
of eigenvalues and eigenfunctions in random Hamiltonian systems are the central issues in 
the applications of RMT. Apart from the nuclear and nuclei fields, it was also employed 
to study the critical statistics of disordered systems with various complex interaction.^"^ 
Recently, resonant properties of composite materials have been studied extensively due to 
the large linear and nonlinear optical responses. In a dielectric network, there exist a 
lot of geometric resonances randomly distributing in the resonant area.^'^^^^ For a specific 
sample, its resonance spectrum is very sensitive to the microstructure. While, for a large 
number of samples given a parameter, i.e., concentration p, the distribution of resonances is 
stable, which implies some general features of resonance spectrum. The aim of this paper is 
to study the level spacing statistics of geometric resonances. 

In this work, a binary network is considered, where the impurity bonds with admittance 
ei are employed to replace the bonds in an otherwise homogeneous network of identical ad- 
mittance 62. The admittance of each bond is generally complex and frequency-dependent. 
All the impurity bonds construct the clusters subspace. For a binary composite, the admit- 
tance ratio h{— ^) of two components has a branch cut on the negative axis when resonance 
happens. Based on the Green's-function formalism (GFF), the eigenvalues(s = and 
s G [0,1]) of Green's-matrix M are solved, the sequence of which forms the resonance 
spectrum. ^^'^^ M maps the geometric configuration of the clusters subspace. Its solutions 
summarize the geometric resonances of the network subject to the external sources and in 
the quasistatic limit. The element of M is defined as Mx,y = Z^zeC(y)(G'x,y — Cx.z), where 
z e C{y) means that the jointing points z and y belong to the impurity metallic cluster 
and are the nearest neighbors. More clearly, M^^y describes the interaction between x and 
y and is closely related to the "environment" or the nearest neighbors of y. In a dilute 
network, there exists a large "environment" difference between two jointing points, hence 
the elements of M distribute more uncorrelatively rather than with some symmetry. How- 
ever, for a percolating network, due to the self-similar structure, x and y have the similar 
"environment" and Mx,y is approximately equal to My^x- It is analogous to the Gaussian 
orthogonal ensemble (GOE), in which the elements of the Hamiltonian matrix must sat- 
isfy Hm,n = Hn,m- So, for a disordered composite, it is expected that the correlation and 
fluctuation of eigenvalues of M have the general features. 

In the following, statistics of resonance level spacing in the disordered binary composites 
is studied intensively on the unfolded scale. For one sample, there are more than 700 levels 
numerically solved by its Green's-matrix M. In the unfolding procedure, we use a fit of the 
third order polynomial to the data. For an arbitrary p, 1000 samples, with totally more 
than 700, 000 levels, are computed for the sake of the ensemble averaging. For a definite 
p within the interval [0.2,0.7], our numerical calculations indicate that the unfolded level 
spacing distribution P(t) and level number variance have the general features. It is 

also shown that the short-range fluctuation P{t) and long-range spectral correlation S^(i^) 
he between the profiles of the Poisson ensemble and GOE. 
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This paper is organized as follows. In the next two sections, for various p, level spacing 
distributions P(t)'s and level number variances E^(L)'s are calculated on the unfolded scale. 
Then, in Section IV, at Pc, crossover behavior of P{t) and S^(L) is obtained, giving the 
finite size scaling of mean level spacing S and mean level number n. Finally, we summarize 
the main results in Section V. 



In order to remove the system-specific mean level density or normalize the resonance level 
spacing of different samples, unfolding procedure is necessary. For a sample, the cumulative 
spectral function C{s) of its resonance spectrum is defined as^ 



where f{s') = Yln=i ^i^' ~ ^n) is the spectral function of levels. When oo, f{s') should 
be smooth. C{s) is the staircase function and is used to count the number of levels. To 
satisfy that f{s') is smooth, s is rescaled by ^ as 



where { is the smooth part of C{s) and Cfi{s) is the fluctuating part of C{s). Fig. 1 
shows the small section of the measured spectrum. In the following, we use the third order 
polynomial to fit the data. However, when p — 0.1 or 0.2, the ninth or eleventh polynomial 
is not high enough to fit the curve because of the large fluctuation of level spacing. So at 
the dilute systems, numerical calculations are not accurate. 

The level spacing distribution P{t) is the probability density for neighboring levels and 
having the spacing t. It is used to describe the short-range spectral fluctuations. On the 
unfolded scale, the function P{t) and its first moment are normalized to unity, /q°° P{t)dt — 1 
and /q°° tP{t)dt = 1. For the uncorrelated or Poisson ensemble Pp{t) = exp{—t), while for the 
GOE, or Wigner-Dyson ensemble, Pwoit) = ^^xp{—^). Fig. 2 displays the level spacing 
distributions P(t)'s for various p within the interval [0.1, 0.7]. In this figure, the distributions 
Pp{t) and Pwoit) are drawn by the dashed and solid hnes respectively. The solid curves 
with the filled or opaque circles represent the critical distributions Prit) at the percolation 
threshold Pc- The inset of Fig. 2 is used to describe the level spacing distributions -P(i)'s 
when p > Pc- In this case, the calculations of the functions P(t)'s are not very accurate 
because of the degeneracy of resonances. For a specific sample, it is impossible to estimate 
where the next level is due to the complex microstructure. However, for the disordered 
composites with a definite p, when the ith level is measured, the spacing between the ith 
and (z + l)th levels satisfies the ensemble averaged distribution P{t), rather than the Poisson 
distribution Pp{t) or GOE case Pwoit). It is seen that the distributions P(t)'s lie between the 
profiles of the Poisson ensemble and GOE for p G [0.1, 0.7]. As discussed in the introductory 
part, the critical Prit) approaches the Pwoit) due to the strong interaction among the 
metaUic bonds. When p is very small, P(t) is close to Pp(i) due to the weak interaction. At 
last, for an arbitrary p, Pit) lies between the functions Ppit) and Prit)- In this figure, it is 
also shown the crossover behavior of Pit). 



II. UNFOLDED LEVEL SPACING DISTRIBUTION 




(1) 



Cis)^C + Cfiis), 



(2) 
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In Fig. 2, wc find that the probabihty of small spacing is much less than the uncorrelated 
Poisson distribution. It means that some repulsion exists between neighboring levels. The 
repulsion in Gaussion ensembles comes from the symmetry of matrix elements. Here, the 
property of matrix elements is described by the Green 's-matrix M of the GFF. So the 
repulsion of resonance levels is caused by the interaction of metallic bonds. The repulsion 
becomes stronger with decreasing \p — pc\. For the GOE, we notice the importance of the 
repulsion law P{t) oc t'^ with f3 = 1 for small spacing.^ While for the composite materials 
with various p, when t 0.0, P(t) > t^. When t > 2.002, here 2.002 is the second 
intersection point of Pp{t) and Pwoit), the long tails of all functions approach the values of 
the Poisson ensemble. 

III. UNFOLDED LEVEL NUMBER VARIANCE 

The above nearest neighbor level spacing distribution contains the information of short 
scales about the resonance spectrum. Long-range correlation is measured by the level number 
variance S^(L), given by^ 

e2(l) c'm e.) > - < cm e.) (3) 

where C{L, ^s) counts the number of levels in the interval [^s, + L] on the unfolded scale. 
The angular bracket denotes the average over the starting points ^s- By the unfolding, 
< C{L, C,s) > should be equal to L. Thus, in the interval of length L, one expects to find the 
L ± ^Ti'^^L) levels on average. For the Poisson spectrum without correlation, one obtains 
Sp(L) = L, while for the GOE, the analytical result is 

= ^(/n(27rL) +^ + 1-^), (4) 

where 7 = 0.5772... is the Euler's constant. Fig. 3 shows the numerical results of level 
number variances S^(L)'s for p e [0.1, 0.7]. The same lines and symbols are used as those in 
Fig. 2, namely, the profiles of the Poisson ensemble and GOE are represented by the dashed 
and solid lines, and the critical level number variance Et^(I/) is plotted by the solid line 
with the circles. It is obvious that the correlation among levels is greater than the Poisson 
case and less than the GOE case. We could not collect the data for p > 0.7 due to the 
degeneracy of eigenvalues. When p = 0.1 or p = 0.2, the level number variance E^(L) is out 
of the boundary of the Poisson ensemble. The reason is that the data can not be fitted very 
well by the third or higher order polynomial in the unfolding procedure. For p e [0.3, 0.7], 
E^(L)'s lie between the profiles Ep(L) and E|,(L). As shown in Fig. 3, the curves almost 
overlap for p = 0.4 and p = 0.6, as well as for p = 0.3 and p = 0.7. So at the percolation 
threshold, the crossover behavior of the level number variance is numerically obtained. 

IV. CROSSOVER AND FINITE SIZE SCALING AT PERCOLATION 

THRESHOLD 

For a percolating network, level spacing distribution limited in the interval s G [0.25, 0.75] 
has been investigated by Luck et. al..^^ Recent numerical calculations indicate that the 
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resonances within [0, 0.25] and [0.75, 1] are important because the high values of the inverse 
participation ratios (IPR) in those regions correspond to the strongly enhanced optical 
responses. In order to study the criticality of level spacing distribution P(t) and level 
number variance S^(L) in a two dimensional network, the duality of level spacing for binary 
model is discussed. We consider an infinite binary network with concentration p. The 
admittance of impurity bonds and matrix bonds is set to ei and 62 respectively. We get 
the admittance ratio h = exje^ and s = The resonance spectrum is given by the 

set {si, S2, ■■■Sn}. Then the set {t} of level spacing can be written as {ti, ^2, •••^n-i} with 
tn-i — Sn — Sn-i- For a large network, G is a typical configuration of the concentration p. 
Binary model is invariant under the simultaneous interchange p <-> 1 — p and ei <-> €2- So 
G' is also a typical configuration of the concentration 1 — p and we get the new admittance 
ratio h' = and s' = 1 — s. The spectrum of resonance is replaced by the set {si', S2') •••Sn'} 
with si' = 1 — si, S2' = 1 — S2 and Sn = 1 — Sn- The new level spacing set {t'} is just the 
original set {t}. So the duality of level spacing exists for the binary model and it is the 
self-dual p — Pc — 0.5. The functions P{t) and are computed based on the same 

set {t} for p and 1 — p. 

In the calculations, we can not find the strict duality of level spacing. One reason is that 
the binary sample is not large enough. Another is that the degeneracy of resonances affects 
the accuracy of P{t) and S^(L) when p > p^. Numerical results indicate the crossover of 
level spacing distribution P{t) and level number variance Y?{L), as shown in Fig. 4 and 
Fig. 5. To characterize this behavior, for each distribution P{t), we compute the value 

^ f^'"'r?^C'^°?y] with to = 0.4729, the first interaction point of PJt) and Pwoit)?'^ In 

Jo y^pW-i^WD\t)\at ' 

this way, r] varies from l[P{t) = Pp{t)\ to 0[P(t) = Pwoit)]. We note that r]T = 0.2162 at 
the transition is close to the value rj^ = 0.215, which corresponds to PA{t) the Anderson 
transition. As observed in Fig. 4, ri{p) is closer to the value 1 of the Poisson ensemble with 
increasing \p — pc\- In Fig. 3, we have found the linear relation of S^(L) with respect to L as 
E^(L) = X* L, where x is called the spectral compressibility.^ When the unfolded number 
L is larger, the value x is 1 for the Poisson ensemble, while for the GOE or Wigner-Dyson 
ensemble, x is approximately equal to 0. For various p, the slopes x's are plotted in Fig. 5. 
At Pc, X = 0.395 is close to 0, rather than the value 1 of the Poisson ensemble. In Fig. 4 
and Fig. 5, we observe the crossover behavior of P{t) and S^(L), as well as the duality of 
level spacing. 

Finite size scaling of mean level spacing 8 and mean level number n are computed when 
the percolating sample is in size from 16 x 16 to 32 x 32. For each case, more than 700, 000 
levels are calculated for the ensemble averaging. The results are shown in Fig. 6 and Fig. 7. 
Finite size scaling laws, 5 = 1.032L~^'^^^ and n = 0.911L^'^^*^, are obtained. Note that here 
the meaning of L is different from that in Fig. 3. The scaling exponents 1.952 and 1.970 
are universal and closely related to the spatial dimension {D — 2) of the network. 



V. CONCLUSIONS 

We have investigated the fiuctuation and correlation of geometric resonance level spacing 
in the random binary composites by RMT. The main conclusions include: 
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1. For a definite the unfolded level spacing P{t) and level number variance S^(-L) have 
the general features. 

2. For an arbitrary p, the short-range spectral fiuctuation P{t) and long-range spectral 
correlation S^(L) lie between the profiles of the Poisson ensemble and critical ensemble. 

3. The functions -Pt(^) and Sf.(L) at the transition approach the profiles of the GOE, 
rather than the profiles of the Poisson ensemble. 

4. The crossover behavior and duality of P{t) and S^(L) are found when p approaches 

Pc- 

5. At Pc, finite size scaling laws, 5 = 1.032L~^'^^^ and n = 0.911I/^'^™, are obtained. 

Statistics of eigenvalues of the Green's-matrix M has been studied and the general sta- 
tistical distributions have been obtained. Statistics of eigenvectors of M, which are closely 
related to the local electric fields, will be published elsewhere. 
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FIGURES 



FIG. 1. The typical cumulative spectral function C(s) of resonance levels. The small part of 
the measured spectrum is shown as the staircase function. The smooth part ^(s) is the third order 
polynomial whose coefficients are found by a fit of the whole measured spectrum. 

FIG. 2. Level spacing distribution P{t) of resonance spectrum on the unfolded scale. Here p is 
ranged at the interval [0.1,0.7]. Dashed and solid lines show -P(t)'s of the Poisson ensemble and 
GOE. The solid lines with filled or opaque symbols show the profiles of various p. 

FIG. 3. Level number variance of resonance spectrum on the unfolded scale. Here p is 

ranged at the interval [0.2,0.7]. Dashed and solid lines show S^(L)'s of the Poisson ensemble and 
GOE. The solid lines with filled or opaque symbols show the profiles of various p. 

FIG. 4. Dependence of rj on the concentration p. The lines are guides to the eyes. 

FIG. 5. Dependence of the spectral compressibility % on the concentration p. The lines are 
guides to the eyes. 

FIG. 6. Finite size scaling of the mean level spacing 6 at Pc- The sample is in size from 16 x 16 
to 32 X 32. 

FIG. 7. Finite size scaling of the mean level number n at Pc- The sample is in size from 16 x 16 
to 32 X 32. 
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L=16,p=0.5, cumulative spectral function C(s). 




r|(p), p is ranged from 0.1 to 0.7. 




The slope D of I^^(L), p is ranged from 0.1 to 0.7. 




Finite size scaling of tine mean level spacing 5, 
the sample size is changed from 16*16 to 32*32. 
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